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We  apply  ultrafast  x-ray  diffraction  with  femtosecond  temporal  resolution  to  monitor  the  lattice 
dynamics  in  a  thin  film  of  multiferroic  BiFe03  after  above-band-gap  photoexcitation.  The  sound-velocity 
limited  evolution  of  the  observed  lattice  strains  indicates  a  quasi-instantaneous  photoinduced  stress  which 
decays  on  a  nanosecond  time  scale.  This  stress  exhibits  an  inhomogeneous  spatial  profile  evidenced  by  the 
broadening  of  the  Bragg  peak.  These  new  data  require  substantial  modification  of  existing  models  of 
photogenerated  stresses  in  BiFe03:  the  relevant  excited  charge  earners  must  remain  localized  to  be 
consistent  with  the  data. 
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Multiferroics  have  a  great  potential  for  application  due  to 
their  possible  coupling  of  ferroelectricity  and  magnetism 
[1-3].  BiFe03  (BFO)  is  one  of  the  few  room  temperature 
multiferroics  today  [4—8],  and  of  these,  the  only  one  that  is 
a  stable  phase.  Its  relatively  small  band  gap  of  approx¬ 
imately  2.7  eV  [9]  renders  BFO  an  ideal  candidate  for 
applications  in  spintronics  and  memory  devices  [5]  with  a 
perspective  for  ultrafast  optical  switching  similar  to  purely 
ferroelectric  [10]  or  magnetic  materials  [11],  The  photo¬ 
voltaic  effect  in  this  complex  material  and  the  underlying 
ultrafast  carrier  dynamics  after  above-band-gap  femtosec¬ 
ond  (fs)  optical  excitation  have  been  studied  thoroughly 
[12-14].  The  photoinduced  currents  in  BFO  lead  to  THz 
emission  [15,16]  and  to  a  photostrictive  response  [17].  All- 
optical  experiments  showed  that  the  rapid  photoinduced 
mechanical  stress  excites  coherent  phonons  [18,19].  The 
dynamics  of  photoinduced  strains  were  directly  and  quan¬ 
titatively  measured  in  a  recent  synchrotron-based  ultrafast 
x-ray  diffraction  (UXRD)  study  with  a  temporal  resolution 
of  lOOps  [20].  Combined  optical  measurements  revealed  a 
linear  dependence  of  the  transient  strain  and  the  number  of 
excited  earners  over  several  nanoseconds  (ns).  This  led  to 
the  conclusion  that  depolarization  field  screening  (DFS) 
including  macroscopic  transport  of  the  carriers  to  the 
surface  and  interface  could  be  the  dominant  stress  gen¬ 
erating  process,  although  the  effect  of  excited  antibonding 
orbitals  was  not  ruled  out  [20]. 


In  this  Letter,  we  report  complementing  UXRD  experi¬ 
ments  at  a  laser-driven  plasma  x-ray  source  (PXS)  in  order 
to  monitor  the  coherent  and  incoherent  lattice  dynamics  in  a 
BFO  thin  film  sample  with  subpicosecond  (ps)  temporal 
resolution  after  above-band-gap  excitation.  We  observe  a 
sound- velocity  limited  evolution  of  the  structural  response 
within  lOps  indicating  a  quasi-instantaneous  stress.  The 
substantial  Bragg  peak  broadening  is  a  direct  evidence  of 
an  inhomogeneous  spatial  stress  profile.  It  appears  quasi- 
instantaneously  and  decays  on  nanosecond  time  scales  as 
reconfirmed  by  new  synchrotron-based  UXRD  data 
recorded  at  the  Advanced  Photon  Source  (APS).  We  obtain 
quantitative  agreement  of  the  transient  peak  shift  and 
broadening  measured  with  both  setups  and  can  firmly 
conclude  that  the  photogenerated  stress  driving  the  film 
expansion  has  a  strongly  inhomogeneous  spatial  profile  in 
the  35  nm  thick  film.  Neither  the  build  up  at  sub-ps  delays 
[21]  nor  the  remaining  spatial  profile  at  ns  delays  of  the 
photoinduced  stress  profile  in  BFO  can  be  determined  by 
all-optical  methods  [22]  alone  and,  hence,  call  for  the 
combination  of  the  applied  UXRD  techniques. 

We  propose  a  model  of  local  charge  earner  displacement 
within  the  BFO  unit  cells  after  above-band-gap  excitation 
resulting  in  a  lattice  distortion  possibly  due  to  the  inverse 
piezoelectric  effect  which  drives  the  expansion  instanta¬ 
neously.  Subsequent  fast  trapping  of  the  excited  charge 
carriers  in  the  film  maintains  the  stress  according  to  the 
optical  excitation  profile  over  several  ns  until  they  decay 
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radiatively  [12].  Our  experimental  study  provides  an 
important  benchmark  for  simulations  of  the  photovoltaic 
response  of  ferroelectric  oxide  materials  [23,24],  which  will 
have  to  predict  strongly  inhomogeneous,  ultrafast,  and  long 
lived  charge  carriers. 

We  investigate  the  very  same  sample  which  was  studied 
in  Ref.  [20].  The  sample  is  composed  of  a  d  —  35  nm  thick 
pseudocubic  (001)  BFO  film  epitaxially  grown  on  a  (001) 
SrTiO^  (STO)  substrate.  The  ferroelectric  polarization 
points  along  the  [111]  pseudocubic  direction  of  BFO 
and  exhibits  a  fourfold  symmetry  with  most  of  the 
polarization  pointing  towards  the  surface  [25].  The  direct 
band  gap  of  this  sample  has  been  determined  to  2.6  eV  [20]. 
At  the  excitation  wavelength  of  A  =  400  nm  the  optical 
penetration  £  =  32  nm  [26]  determines  the  excitation 
profile  following  the  Lambert-Beer  law. 

The  UXRD  setups  at  the  PXS  and  the  APS  have  been 
described  elsewhere  [27-29].  The  PXS  provides  a  temporal 
resolution  below  200  fs  at  an  x-ray  photon  energy  of 
8.047  keV  (Cu  Ka)  and  is  operated  in  a  convergence- 
correction  mode  [30].  The  x-ray  and  UV  footprints  on  the 
sample  have  diameters  of  approximately  300  //m  and  1  mm 
(FWHM),  respectively.  The  UV  pump  beam  is  p  polarized 
and  incidents  at  40°  from  the  surface.  The  synchrotron 
based  setup  provides  much  higher  stability  for  long  term 
measurements  in  the  ns  range,  while  excitation  and  probing 
conditions  are  very  similar  [20]. 

Figure  1(a)  shows  the  static  rocking  curve  of  the  002 
pseudocubic  Bragg  peaks  of  the  BFO  and  STO  substrate  as 
measured  by  the  PXS  together  with  a  transient  rocking 
curve  at  t  =  10  ps  delay.  The  high  crystalline  quality  of  the 
film  is  evidenced  by  the  the  static  rocking  curve  which 
coincides  with  the  dynamical  x-ray  simulation  for  a 
d  =  35  nm  thick  perfect  BFO  film  on  STO  including  the 
instrumental  resolution.  Figure  1(b)  shows  the  transient 
rocking  curves  for  the  early  delays  from  —10  to  25  ps 
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FIG.  1  (color  online),  (a)  Rocking  curves  of  the  002  pseudo¬ 
cubic  Bragg  reflections  from  the  BFO  layer  and  STO  substrate 
measured  by  the  PXS.  The  solid  black  line  is  a  simulation  of  the 
static  rocking  curve,  (b)  Measured  transient  rocking  curves 
(diffracted  intensity  in  logarithmic  scale).  The  green  dashed  line 
indicates  the  extracted  center  of  the  BFO  peak. 


with  the  photoexcitation  occurring  at  t  =  0.  The  peak  shift 
of  the  BFO  002  pseudocubic  reflection  measures  the 
average  out-of-plane  strain  (z  direction)  in  the  layer: 
—A qz(t)  <x  (e(z,  t))  .  The  observed  shift  to  smaller  qz 
corresponds  to  an  ultrafast  expansion  Ac/ c  «  0.5%  of 
the  BFO  film  along  the  surface  normal  without  any 
contraction  features  which  were  observed  for  the  ferro¬ 
electric  material  PbTi03  (PTO)  [31]. 

Figure  2  shows  the  linear  dependence  of  the  transient  peak 
shift  on  the  absorbed  fluence  for  selected  delays  ranging 
from  t  =  0.01  to  15  ns  and  confirms  that  data  from  the  PXS 
and  APS  setups  quantitatively  agree  within  a  reasonable 
40%  recalibration  of  fluences  between  the  two  laboratories. 
The  linear  fluence  dependence  suggests  that  the  origin  of 
stress  is  the  same  for  early  (ps)  and  late  (ns)  delays. 

The  transient  BFO  peak  shift  and  width  are  plotted  in 
Fig.  3.  The  smaller  peak  shift  of  the  APS  data  within  the 
first  50  ps  originates  from  the  limited  temporal  resolution 
(lOOps)  of  the  synchrotron-based  experiments.  The  decay 
of  the  shift  is  fitted  by  an  exponential  function  with  a  time 
constant  ~  2.29  ±  0.14  ns  including  an  offset  for 
very  late  time  scales. 

The  width  w  inversely  depends  on  the  number  of 
scattering  atomic  layers  (size  broadening)  and  on  the 
inhomogeneous  strain  fields  within  the  film  (strain  broad¬ 
ening)  [32,33].  We  can  neglect  mosaic  broadening  for  the 
high-quality  BFO  film  [34]  and  for  the  synchrotron-based 
APS  experiments  instrumental  broadening  is  negligible,  too. 

Although  the  link  between  the  peak  profile  and  the 
spatial  strain  profile  is  complicated,  we  can  assume  that  the 
change  of  the  width  Aw  depends  linearly  on  the  spatial 
variation  A Ze{z,t)  of  the  strain:  Aw(f)  oc  A ,e(z,  f).  If  the 
layer  is  homogeneously  strained  (Aze  =  0),  no  additional 
peak  broadening  is  observed  (Aw  =  0).  The  transient  width 
Aw  reveals  a  significant  inhomogeneous  strain  profile  in 
the  BFO  layer  over  the  whole  observed  time  scale.  The 
exponential  fit  in  Fig.  3(b)  results  in  comparable  decay 
times  for  the  width  and  shift:  « 2.31  ±  0.92  ns. 

This  directly  indicates  that  the  spatial  strain  profile 
does  not  equilibrate  within  the  thermal  relaxation  time 


0.5 

0.4 

0.3 

0.2 

0.1 

0.0 


FIG.  2.  Shift  A qz  of  the  002  pseudocubic  BFO  Bragg  peak  and 
corresponding  average  strain  (e(r))_  measured  as  a  function  of 
fluence  at  different  delays.  The  dashed  and  dotted  lines  represent 
linear  fits. 
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FIG.  3.  The  transient  shift  A qz  and  corresponding  average 
strain  (e(t))z  (a)  as  well  as  the  relative  change  of  the  width 
Aw/w0  (b)  of  the  002  pseudocubic  BFO  Bragg  peak.  The  APS 
data  were  measured  at  an  absorbed  fluence  of  F  =  2.4  mJ/cm2 
and  the  PXS  data  was  scaled  according  to  the  fluence  calibration. 
The  black  solid  line  represents  an  exponential  fit  of  the  APS  data. 

rth  —  d2pc/K  —  850  ps  of  the  d  =  35  nm  thick  BFO 
layer  [35]. 

For  delays  t  20  ps  larger  than  the  time  it  takes  strain 
waves  to  travel  twice  through  the  thin  film  at  the  speed  of 
longitudinal  sound  nBFO  =  3.5  nm/ps  [19],  the  observed 
shift  A qz  is  not  only  proportional  to  the  average  strain  but 
also  to  the  transient  stress  according  Hook’s  law 
—A q.  oc  (e(z,  t))  <x  (a(z,  t))v  The  similarity  of  the  decay 
times  for  the  peak  shift  and  broadening,  moreover,  suggests 
that  the  transient  stress  o(z,  t)  can  be  approximated  by  a 
time-invariant  spatial  stress  profile  f(z)  that  decays  in 
amplitude  A(t);  i.e.,  no  transport  processes  are  relevant  to 
the  driving  stress:  a(z,t)  =  A(t)  x  f(z).  We  can  then 
directly  link  the  peak  shift  A q:(t)  to  the  amplitude  A(t) 
by  —A qz(t)  oc  (a(z,t))z  oc  A(t)  x  (f(z))z,  as  well  as  the 
peak  width  by  Aw(f)  oc  A{t)  x  A Zf(z). 

To  elucidate  the  origin  of  the  photoinduced  stress  at  very 
early  delays  (f  <  20  ps),  where  the  stresses  cannot  be 
calculated  via  Hook’s  law,  we  simulate  the  BFO  peak  shift 
from  a  ID  lattice  dynamics  simulation  [36,37]  of  the  strain 
profile  e(z,  t)  for  the  given  stress 

<j(z,  t)  oc  A(t)  x  f(z)  =  H(t)(  1  —  ye~,lT^)  x  e~z^.  (1) 

We  explicitly  use  the  optical  penetration  depth  £  =  32  nm 
and  assume  a  time-dependent  rise  A{t)  of  the  stress  which 
includes  a  quasi-instantaneous  stress  approximated  by  the 
Heaviside  step  function  H(t)  and  an  additional  stress 
component  growing  with  the  time  constant  rrise.  The 
transient  stress  is  plotted  in  Fig.  4(a)  for  several  relative 
strengths  y  of  the  instantaneous  and  delayed  stress  com¬ 
ponents  [38].  We  apply  dynamical  x-ray  theory  to  calculate 
the  according  rocking  curves  from  the  simulated  strain 
e(z,  t)  in  order  to  extract  the  transient  peak  shift  in  the  same 
way  as  for  the  experimental  data  [37].  We  clearly  obtain  the 
best  fit  to  the  experimental  peak  shift  for  y  —  0.  In  the 


FIG.  4.  The  black  circles  in  (b)  show  the  transient  shift  of  the 
002  pseudocubic  BFO  Bragg  peak  extracted  from  the  data  shown 
in  Fig.  1(b).  The  black  solid  line  represent  the  simulated  peak 
shift  for  an  instantaneous  stress.  The  simulated  peak  shift  for  a 
semi-instantaneous  stress  (gray  dashed  line)  and  a  fully  diffusive- 
like  stress  (gray  dotted  line)  both  with  a  time  constant 
TriSe  =  2  ps,  see  text.  The  corresponding  stress  amplitudes  A(t) 
are  plotted  in  panel  (a).  The  gray  Gaussian  at  t  =  0  in  panel  (b) 
indicates  the  temporal  resolution  of  200  fs  of  the  PXS  setup. 

examples  shown  in  Fig.  4,  we  used  the  time  scale 
rrise  —  2  ps;  however,  we  have  cross-checked  this  state¬ 
ment  with  additional  simulations  using  longer  time  scales 
rnse  f°r  the  delayed  stress  and,  also,  for  different  spatial 
stress  profiles.  In  essence,  this  proves  that  the  dominant 
contribution  to  the  stress  is  instantaneous  and  spatially 
inhomogeneous. 

We  want  to  emphasize  the  importance  of  several  key 
parameters  entering  the  simulation:  (i)  The  temporal  over¬ 
lap  of  x-ray  probe  and  UV  pump  pulses  ( t  =  0)  was 
determined  independently  with  an  accuracy  of  approx. 
±100  fs  [21].  (ii)  The  film  thickness  was  determined 
experimentally  by  XRD  measurements  to  d  =  35  nm. 
(iii)  We  determined  the  longitudinal  acoustic  sound  veloc¬ 
ity  ^bfo  =  3.5  nm/ps  independently  by  an  oscillation  of 
the  layer  thickness  after  photoexcitation  [38]  (not  shown) 
with  good  agreement  to  literature  values  [19]. 

The  DFS  model  proposed  as  origin  of  the  fast  lattice 
expansion  in  BFO  [20]  and  in  ferroelectric  PTO  [31] 
requires  free  charge  carriers  to  travel  from  within  the  bulk 
material  to  the  interfaces  of  the  layer  to  screen  the 
depolarization  fields.  The  ferroelectric  polarization 
increases  while  the  carriers  drift  with  velocity  vd.  The 
inverse  piezoelectric  effect  would,  thus,  yield  a  finite  stress 
rise  time  xc  which  is  required  for  the  earners  to  propagate 
across  the  film.  The  typical  charge  carrier  mobility  in 
ferroelectrics  is  between  //  =  0. 1-3.0  cm2/Vs  [39]  and 
typical  internal  electric  fields  are  in  the  range  of  E « 
200  kV/cm  [40].  Even  for  a  high  charge  carrier  mobility  of 
f,i  <  3  cm2/Vs,  this  leads  to  rc  =  d/vd  =  d/ iaE  >  5  ps  for 
the  BFO  sample  which  contradicts  the  quasi-instantaneous 
stress  that  is  required  to  drive  the  ultrafast  lattice  dynamics, 
cf.  Fig.  4.  We  would  further  expect  a  spatially 
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homogeneous  stress  profile  from  the  DFS  model  due  to  the 
capacitorlike  geometry  in  the  film,  which  is  incompatible 
with  the  significant  long-lived  peak  broadening  shown  in 
Fig.  3(b).  We  conclude  that  DFS  due  to  diffusion  of  charge 
carriers  across  the  whole  film  thickness  cannot  be  a 
dominant  process  in  BFO. 

Unfortunately,  the  fluence  dependent  study  (Fig.  2) 
cannot  differentiate  between  alternative  stress  mechanisms, 
since  all  of  them  essentially  depend  linearly  on  the  fluence. 
In  the  DFS  model,  this  is  only  true  for  low  fluence,  since 
the  process  saturates  when  the  depolarization  field  is  fully 
screened  [41,42].  For  PTO  [31],  the  saturation  occurs  at 
excitation  fluences  of  approximately  1  mJ/cm2.  In  our 
experiment,  the  strain  in  BFO  remains  linear  up  to  absorbed 
fluences  of  more  than  2  mJ/ cm2  (Fig.  2)  with  comparable 
strain  levels.  At  even  higher  fluences.  the  BFO  film  starts  to 
degrade.  Since  only  a  very  small  number  of  the  excited 
carriers  is  sufficient  to  drive  the  stress  in  the  DFS  model 
[20],  the  absence  of  strain  saturation  effect  disfavors  the 
DFS  model  as  well,  possibly  due  to  the  presence  of  a  skin 
layer  at  the  BFO  surface,  which  might  change  the  ferro¬ 
electric  boundary  conditions  [43]. 

The  thermal  contribution  to  the  stress  can  be  quantitatively 
estimated:  The  electrons  in  the  BFO  film  exhibit  an  excess 
energy  of  approximately  0.5  eV  in  the  conduction  band  when 
they  are  excited  with  3.1  eV  photons.  This  excess  energy  is 
transferred  rapidly  to  the  lattice  via  electron-phonon  coupling 
within  1  ps  [12,20]  resulting  in  a  fast  temperature  increase  of 
the  lattice.  We  calculate  this  temperature  jump  as  AT  = 
44  K  averaged  over  the  BFO  film  thickness  for  an  absorbed 
fluence  of  F  —  2.4  mJ / cm2  taking  into  account  the  internal 
refraction  of  the  pump  light  in  the  BFO  layer.  This  corre¬ 
sponds  to  a  maximum  thermal  strain  of  e  —  44  K  x  1.84  x 
10~5  K  1  =  0.08%  which  accounts  only  for  a  small  fraction 
of  the  peak  shift  [Fig.  3(a)].  A  numerical  simulation  [37]  of 
the  heat  diffusion  assuming  bulk  values  [20]  yields  the  time 
scale  of  5  ns  for  reducing  the  average  BFO  temperature  rise 
to  6  K.  Even  the  artificial  introduction  of  a  very  large  thermal 
interface  resistance  [44]  keeping  the  heat  in  the  BFO  layer 
would  only  explain  about  half  of  the  observed  shift 
for  t  >  10  ns. 

The  dominant  mechanism  right  after  excitation  leads  to 
an  immediate  and  inhomogeneous  lattice  expansion,  and  is 
still  important  after  10ns.  The  optical  excitation  transfers 
electrons  into  initially  unoccupied  orbitals  above  the  band 
gap  where  they  can  contribute  to  the  lattice  dynamics  by 
deformation  potential  interaction  [45].  In  contrast  to  our 
experimental  observation,  this  effect  would  lead  to  a 
contraction  of  the  BFO  film  referring  to  Eq.  (17)  of 
Ref.  [22]  and  taking  the  pressure-dependence  of  the  band 
gap  in  BFO  into  account  [46,47].  However,  at  the  same 
time,  the  redistribution  of  electron  density  within  the  BFO 
unit  cell,  also  referred  to  as  a  shift  current  [24],  influences 
the  local  ferroelectric  polarization,  as  evidenced  by  THz 
emission  [15,16].  The  change  in  amplitude  and  sign  of  the 
local  ferroelectric  polarization  depends  on  the  involved 


orbitals  and,  thus,  also  on  the  optical  excitation  wavelength 
as  well  as  on  the  specific  sample  configuration  [23,24].  For 
our  case,  the  photoinduced  change  of  the  local  ferroelectric 
polarization  induces  a  stress  via  the  inverse  piezoelectric 
effect  which  drives  the  ultrafast  lattice  expansion. 
Compressive  stress  due  to  deformation  potential  interaction 
is  overcompensated  for  the  entire  excited  state  lifetime  of 
the  charge  carriers. 

Theoretical  and  experimental  studies  provide  evidence 
for  the  optical  generation  of  p-d  charge  transfer  (CT) 
excitonic  states  in  BFO  [12,48].  These  states  modify  the 
local  ferroelectric  polarization  and  also  couple  strongly  to 
the  structure  of  the  BFO  unit  cell  [48].  The  strong 
absorption  edge  smearing  of  the  2.6-2. 8  eV  band  [26,49] 
and  an  additional  weak  2.4  eV  band  superimposed  on  its 
tail,  measured  by  static  spectroscopy,  point  to  a  CT 
instability  as  well  as  self-trapping  of  p-d  CT  excitons 
and  the  nucleation  of  electron-hole  droplets  [48].  Because 
of  the  rapid  trapping  of  the  CT  excitonic  states,  the  initial 
spatial  excitation  profile  is  maintained  and  the  temporal 
dependence  of  the  lattice  strains  is  solely  determined  by 
electronic  recombination  and  not  by  additional  diffusion 
processes.  We  note  that  the  optical  excitation  is  not 
sufficient  to  drive  the  system  into  a  different  structural 
phase.  A  transient  structural  phase,  which  can  arise  from 
electronic  excitation  or  lattice  distortion  [50],  is  not 
observed  in  our  diffraction  study. 

In  conclusion,  we  applied  UXRD  experiments  with  sub- 
ps  temporal  resolution  to  monitor  transient  lattice  dynamics 
in  a  multiferroic  BFO  thin  film  after  above-band-gap 
photoexcitation.  The  peak  shift  reveals  a  rapid  expansion 
that  is  only  limited  by  the  sound  velocity,  indicating  a 
quasi-instantaneous  photoinduced  stress.  The  peak  broad¬ 
ening  indicates  a  strongly  inhomogeneous  spatial  stress 
profile  for  ps  up  to  ns  delays,  excluding  thermal  stresses  as 
the  dominant  process.  We  propose  a  model  of  a  local  charge 
carrier  displacement  within  the  BFO  unit  cell  after  photo¬ 
excitation  leading  to  an  instantaneous  stress  due  to  the 
inverse  piezoelectric  effect.  The  fast  trapping  of  the 
involved  charge  carriers  maintains  the  spatial  excitation 
profile  until  they  decay  radiatively  on  a  ns  time  scale.  We 
believe  that  this  information  is  essential  for  testing  theo¬ 
retical  models  that  can  also  distinguish  the  contribution  of 
the  piezoelectric  effect  and  of  the  deformation  potential. 
The  subtle  distinction  between  these  processes  could,  in 
principle,  be  given  by  more  nontrivial  UXRD  experiments 
where  an  additional  (e.g.,  electronic)  control  over  the 
polarization  is  implemented. 
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